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CN I Abstract 

Oh| Belle has observed surprisingly copious production of ip{3770) in B meson decay 

^ ; B ^ tp{3770)K, of which the rate is comparable to that of S ^ 1/^(3686)7^. We 

■ study this puzzling process in the QCD factorization approach with the effect of S-D 
mixing considered. We find that the soft scattering effects in the spectator interactions 

CN I play an essential role. With a proper parametrization for the higher twist soft end- 

' point singularities associated with kaon, and with the S-D mixing angle 9 = —12°, the 

. calculated decay rates can be close to the data. Implications of these soft spectator 

I effects to other charmonium production in B exclusive decays are also emphasized. 

o : 

\0 ■ PACS numbers: 13.25.Hw; 12.38.Bx; 14.40. Gx 

o ■ 

1 Introduction 

Oh! 

^ ' The V^(3770) is the lowest lying charmonium state above the open charm DD threshold. It is 
^ ■ expected to be predominantly the l^Di charmonium state with a small admixture of the 2^5"! 
K~j ' component. The ^^(3770) is of great interest in recent studies of charmonium physics. There 
rS I are a number of new measurements and related theoretical issues about the '?/'(3770), e.g. 
I the non D — D decays including charmonium transitions and decays to light hadrons^l2lEl 
(see also jl]), the radiative transitions to the P-wave charmonia'S', the S — D mixing, and 
the discussions about the well known pvr puzzle in J/ip and ?/'(3686) decays (see, e.g., jElIZI)- 
In this paper, we will focus on another interesting issue about the iIj{3770). That is the 
■?/'(3770) production in the B meson exclusive decay B —>■ iIj{3770)K, of which the rate is 
found by Belle to be surprisingly largejS], even comparable to that of i? ^ iIj{368Q)K, and 
it might seemingly indicate that this result suggests a large amount of S-D mixing in the 
i'{S770)\^. But this apparently needs a careful examination. 

It is generally believed that if the virtual charmed meson pair components are neglected 
the two states '?/'(3686) and '?/^(3770) can be approximately expressed as 

I V') =1 ^^(3686)) = cose I cc{2^Si)) + sin9 \ cc(l^Di)), 

I =1 ^(3770)) = cosO I cc(l^Di)) - sine \ cc{2^Si)). (1) 

The S-D mixing angle has been estimated by using the ratio of the leptonic decay widths 
[Oj of '?/'(3686) and '?/'(3770). Nonrelativistic potential model calculations give two solutions: 
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9 ^ -10° to -13° or 6* ^ +30° to +26° 0111 CUl- The small mixing angle is compatible with 
the results obtained in models with coupled channel effects fUE] and is favored by the El 
transition ip' —> •yXcj data also^Oj. 

Belle Collaboration |H] has observed ?/'(3770) in the B meson decay — > ■ip{3770)K^ 
with a branching ratio, 

Br(5+^^"ir+) = (0.48 + 0.11 + 0.07) X 10"^ (2) 

which is comparable to that of ■?/^(3686) 0, 

By{B+-^^'K+) = (0.66 + 0.06) X 10"^ (3) 

This is quite surprising, since conventionally the '?/^(3770) and '?/^(3686) are regarded as pre- 
dominantly the l^Di and 2^5"! cc states respectively, and the coupling of l^Di to the cc 
vector current in the weak decay effective hamiltonian is much weaker than that of 2^ Si 
in the naive factorization approach ^3]. If this experimental result is really due to a large 
S-D mixing, as suggested in |H|, then it is found in Ref.jTl] that an unexpectedly large S-D 
mixing angle 6 = +40° would be required by fitting the observed ratio of 5 — >■ '?/'(3770)ii' to 
B iIj{368Q)K decay rates, when the D wave contribution is neglected. This is in serious 
contradiction with all other experimental and theoretical studies, and, in particular, with the 
newly measured El transition rates for iIj{3770) lXcj{J = 0, 1,2), for which the CLEO 
resultspl are 172 + 30, 70 + 17, < 21 KeV respectively for J=0,l,2 whereas the corresponding 
calculations are 386,0.32,66 KeV for 6 = -40° and 52,203,28 KeV for 9 = +40°tl4j. So, 
based on the naive factorization, the use of large S-D mixing to explain the Belle data for 
B — s> iIj{3770)K should be ruled out. The next question is, can we explain the Belle data by 
considering the nonfactorizable contributions to these decay rates? 

In the following, we will study this problem in the QCD factorization approach |15| I16| ITT] 
including nonfactorizable contributions. We will first give the decay rate of -B — ip{3770)K 
based on the assumption that '?/'(3770) is a pure D-wave charmonium state. Then we take 
the S-D mixing into account. Finally, we will consider the higher twist effects. 



2 B ^ ilj{3770)K decay in QCD factorization 

The effective Hamiltonian for this decay mode is written as[TH] 

/-y 10 

Hes = (VcbKliCiOi + C,02) - VttV* Y: aO,) . (4) 
V 2 i=3 

Here Cj's are the Wilson coefficients which can be evaluated by the renormalization group 
approach JH] and the results at fi = 4.4 GeV are listed in Tab. 1. The relevant operators 
Oi in i^efj are given by 

Oi = {Sabi3)v-A ■ {C(sCa)v-A, O2 = {Saba)v-A " (C/3C/3)y_A, 

^3(5) = {Saba)v-A ' Y{Qi3Ql3)v-AiV+A), (^4(6) = {Sabf3)v-A " Y{Qi3Qa)v-AiV+A), (5) 

q 1 
3 3 

^7(9) = ■7:{Saba)v-A " Y<^q{Q/3(ll3)v+A{V-A), Os(10) = -{Sabf3)v-A ' 'Yeq{qf3(la)v+A{V-A)- 
^ a ^ n 



2 





Ci 


C2 


C3 


C4 


C5 


Ce 


LO 
NDR 


1.144 
1.082 


-0.308 
-0.185 


0.014 
0.014 


-0.030 
-0.035 


0.009 
0.009 


-0.038 
-0.041 



Table 1: Leading-order(LO) and Next-to-leading-order(NLO) Wilson coefficients in the 
NDR scheme (See Ref.|IHI) with = 4.4 GeV and & = 225 MeV. 



We treat the charmonium as a color-singlet non-relativistic cc bound state. Let be 
the total 4-momentum of the charmonium and 2g^ be the relative 4-momentum between c 
and c quarks. For D-wave charmonium, because the wave function and its first derivative at 
the origin vanish, 7^i(0)=0, 7^']^(0)=0, which correspond to the zeroth and the first order in 
q, we must expand the amplitude to second order in q. Thus we have (see, e.g.. 
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x5(g°-^)^;M(9)Tr[On0)Pi5.(p,0) + O"(0)P,l(p,0) 

+O%0)P,%ip,0) +O(0)At(p,0)], (6) 

where 0{q) represents the rest of the decay matrix element. The spin-triplet projection 
operators Pis^ (p, q) is constructed in terms of quark and anti-quark spinors as 

/ 3 'p J) 
PisAp^ q) =\ — Y.'"i7^-(l,S2)u{-^ + q,si){si;s2\lS,), (7) 



and 



oqa OqaOqfs 



(J=0; 



PisMO) - l.=o, PisMO)- ^^^^^^ |,=o. (8) 

After g° is integrated out, the integral in Eq.® is proportional to the second derivative 
of the D-wave wave function at the origin by 



,,q''q'r2M = V^7^;(0), (9) 
[zTij-^ V on 

where e'^ is the polarization tensor of an angular momentum-2 system and the value of 
7^£,(0) for charmonia can be found in, e.g., Ref.j2D]- 

The spin projection operators Pi5^ (p, 0) , P^g^p, 0) andP{^(p, 0) can be written as[T^ 



PisMO) = d—f{S.W + M), (10) 



PisSp^O) = J^s[fiS.W+M)r+rfiS.W+M)], (11) 



where we have made use of the non-relativistic approximation for the charmonium mass 
M ~ 2m. Here m is the charmed quark mass. 

As for the hght meson kaon, we describe it relativistically by hght-cone distribution 
amphtudes (LCDAs) [T7] up to twist-3 leveh 



4 Jo 



y<{H,My) - f^Ki^my) - <y,.v'{z, - z.y ^)}^^, (12) 



where y and y = 1 — y are the momentum fractions of the s and d quarks inside the K 
meson respectively. Here the chirally enhanced mass scale = /{^sifJ') + ^difJ^)) 
is comparable to mb, which ensures that the twist-3 spectator interactions are numerically 
large, though they are suppressed by l/rrih. The twist-2 LCDA (pxiy) and the twist-3 ones 
(p^xiy) ^"^^ ^'xiy) symmetric under y ^ y in the limit of SU(3) isospin symmetry. In 
practice, we choose the asymptotic forms for these LCDAs, 

My) = rK{y) = Qy{'^-y), = i- (is) 

In the naive factorization, we neglect the strong interaction corrections and the power 
corrections in Aqcd/'"^6- Then the decay amplitude can be written as 



iMo = -fDm^" {2pB ■ enF,{ml„)^[v,,V:,{C, + ^) - Vt.V.UCs + ^ + ^5 + ^)J , (14) 



where is the number of colors. We do not include the effects of the electroweak penguin 
operators since they are numerically small. The form factors for B ^ K are given as 



,2 ^2 



{K{pK)\-si,h\B{PB)) = [{pb+Pk), - '^^^-p^P,]F,{p') + !!^^^p,Fo(p^), (15) 



where p = Pb —Pk is the momentum of if) with p^ = rn^,, . The kaon mass will be neglected 
in the heavy quark limit and we will use the approximate relation Fq {m'^n ) / Fi {fn'^„ ) = 1 — r 
\21\ I22j. where r = m?^ii/m\, to simplify the amplitude in our calculations. 

As we can easily see in Eq. dH)), this amplitude is unphysical because the Wilson coeffi- 
cients depend on the renormalization scale /x while the decay constant and the form factors 
are independent of /i. This is the well known problem with the naive factorization. How- 
ever, if we include the order corrections, it turns out that the /i dependence of the Wilson 
coefficients is largely cancelled and the overall amplitude is insensitive to the renormaliza- 
tion scale. Taking the nonfactorizable order as strong-interaction corrections in Fig. ^ into 
account, the full decay amplitude for B —>■ K within the QCD factorization approach is 
written as 

Gi 
V2 

where the coefficients {i = 2, 3, 5) in the naive dimension regularization (NDR) scheme 
are given by 

a, = -{C, + ^) + ^^C,{-12ln^ + 2 + f, + fny 

Nc 4:71 Nc ^ [1 ' 



%M = /Dm^"(2pB ■ £*)F,{ml„)^\v,kV:^a2 - Vt^VZia^ + a,)\, (16) 
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(a) 



c c c c 




(e) (/) 

Figure 1: Feynman diagrams for nonfactorizable corrections to B ip" K decay. 



-(C, + &) + ^^ft(-121u!^ + 2 + /, + /„ 



(17) 



as = - 

a5 = -iC. + ^) - f^^C:6(-12 In ^ - 10 + /, + . 
The function fi in Eq.lfTTj) is calculated from the four vertex diagrams (a,b,c,d) in Fig.[TJ 



dx 



dy 



-61n[(x + |)(x + ^)|((r-l)a: + ^)] 



-2r?/(l - ?/)( 
-2( 



x + f)(x + f) f((r-l)x + f) 
(1 + r)x — r(2 — y) (r — \)x — r(2 — ?/) 



(r-l)x + f 



— 2r(l — r)x2/^ 



■(x + f)(x + f)2 |((r-l)x + f)- 



r) 



+ -r(l — r) X y { 



+ 



where r = th? ,, /m\. 



(18) 



The function /// in Eq. (jl7p is calculated from the two spectator interaction diagrams 
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(e,f ) in Fig. ^ and it is given by 



fKf. 



B 



X 



C 



dy(t)K{y) 



10(l-r)(l-|/) (l-2/)2(l_r)^ 



(19) 



where 0b is the hght-cone wave functions for the B meson. The spectator contribution 
depends on the wave function 0^ through the integral 



/■ 

JO 



(20) 



Since (f>B{0 is appreciable only for ^ of order Aqcb/^b, ^b is of order Aqcd- We will follow 
Ref. to choose ~ 300 MeV in the numerical calculation. 

It is easily seen from ()19p that there is logarithmic end point singularity in the integration 
over y when y ^ 1. It breaks down the factorization even at the leading twist level. It 
implicates that the soft mechanisms may be important to this decay mode. To estimate 
these soft effects, we simply parameterize the end point singularity as 



X 



^ ^ - In [— 



(21) 



where A/j ~ 500 MeV is the typical momentum scale associated with the light quark in 
the B meson. Furthermore, since the virtuality of the gluon exchanged between the spec- 
tator quark and the charm (or anti-charm) quark is Ah^b, we should multiply a factor 
as{y/ Afi'mb)Ci{\/ AhiTib) / {as{n)Ci{n)) to /// in Eq (fT7|) . where /i ~ m;, is the scale at which 
we evaluate those vertex corrections. 

The decay constant fn is calculated through the potential models, 



D 



27rm„ 



m 



(22) 



For numerical analysis, we choose Fi(m^„) = 0.97 23J and use the following input parame- 



ters: 



nib 
fB- 



: 4.8 GeV, niB = 5.28 GeV, m^" = 3.77 GeV, 
216 MeV|24j, fx = 160 MeV. 



(23) 



Then we get the branching ratio: Br(i? ip" K) = 1.13 x 10 ^. The theoretical calcula- 
tion is about 40 times lower than the experimental data 



3 B ^ ip'K decay 

The calculation of the branching ratio for B —>■ ip'K decay is similar to that for B — > il)"K. 
If one treats ip' as a pure 25'-state, the only modification needed to do is to expand the decay 
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amplitudes to zeroth order in the g-expansion. Thus the amphtudes will be proportional to 
the S-wave wave function at the origin through the integration 



(24) 



The full decay amplitude for B ^ ip K within the QCD factorization approach is written 



as 



iM 



\ 



R2s{0)m^>{2pB-e*)F,{m: 



G p 

'71 



ycbV*ao - VtbV*(a. + a 



(25) 



where the coefficients a[ {i = 2, 3, 5) in the naive dimension regularization(NDR) scheme are 
given by 



(C5 + 



rrib 



iih 



(26) 



^5(12 In— + 10 + /, + /,, 



Again, the vertex corrections associated with Fj are evaluated at renormalization scale n 
nib and the spectator interactions associated with Fjj are evaluated at ^/ Ahmb. 
The function /, and /,, in Eq. ()2(ij] have following forms, 



+ 4- 



dx 
2x(l 







dy [61n[(x + |)(x + y)|((2; 



l)^ + f)] 



zy — {1 — z)x 



in 



x + f i((z-l)x+f, 

fKfB ,>b(0 



II 



N,mlF,{ml,){l-z) Jo 



d^- 



dy- 



(y) 



(27) 



where z = m?^, /m?^ and Fi{m?^,) = 0.91. One can easily get the functions in ()27|1 by using the 
known results of i? ^ J/ipK in Ref. [21], where J/ip is described by LCDAs. We only need 
to replace the decay constant fj/^ by f2s = . / ^ ^ -^25(0) and choose the non-relativistic 



limit form 4>nr{u) = 6{u — 1/2) for LCDAs of J/ip as in Ref. 

According to Eq. (1) we can write down the ratio of the decay rates directly: 



R 



Br{B i)"K) 
Br{B ij'K) 



-iM. X sinO + iM. x cos9 
iM' X cos6 + iM X sinO 



(28) 



The ratio determined by experimental data is ~ 0.72. Comparing it with our calcula- 
tion, we can find the mixing angle is: 6 = —26° or 6' = +59°. But the absolute branching 
ratio of B —>■ ip K is 5.9 x 10~^, which is still about one order of magnitude lower than 
experimental data in Eq. 
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4 Higher twist effects and end point singularities 



^0 Ah 



In last two sections, we have only considered the leading twist spectator interactions. Gen- 
erally, the contributions arising from higher twist LCDAs of K meson will be suppressed 
by powers of l/rrib. However, as we have mentioned, the chirally enhanced scale hk ~ iTib 
in ()12|1 ensures that the twist-3 contributions may be numerically large. It was discussed 
several years ago that these contributions may play important roles in the process of B me- 
son to S-wave charmonia decays [22]. Here we consider the higher twist effects in D-wave 
charmonium production as well. 

The distribution amplitude of kaon to twist-3 have been given in ()12j) , then we can find 
the twist-3 modifications to /// and fjj to be 

N, mlFi{m^^„) {1 - ry Jo ^ ^ Jo ^ 6 \r - 1) 10 y^^' 

Nc m%Fi{m^,) [1 — z)'' Jo t, Jo Q 

where = 2/ii^/mfe. Here, we can see there exist logarithmic end point singularities in 
both fji and /pj. More seriously, there emerges linear singularity in function ffj and we will 
parameterize it just like what we have done for the logarithmic ones: 

''^ dy TUB 
lo y^ ~ Ah 

It is implicit that these singularities can be regularize by the gluon or light quark off- 
shellness of order when we use ()2H) and ()30p. So when the offshellness is negative, the 
logarithmic singularity will receive large complex contributions from the implicit pole in the 
region of integration. They are common effects in soft rescattering processes. Then following 
|2S], we rewrite (pT|l as 

A-./'*=ln(!^)+t. (31) 

Jo y \Ah/ 

where t is a complex free parameter and we choose varying from 3 to 6 as suggested in 
[22] • Setting |t| = 4.5, < 5 < vr, and the S-D mixing angle 9 = —12°, we can get the 
branching ratio curves of i? ^ iIj{3770)K and B — > iIj{3686)K, which are shown in fig|21 

From fig|21we see that in the region with small 6 the branching ratios are not very sensitive 
to the value of 6. With a value of, say 6 = vr/S, the branching ratios of -B — tp{3770)K and 
B 'iIj{3686)K are found to be: 

Br(5+ ^ ^"K+) = 2.68 x 10"^ 

Br(5+ ^ ^'K^) = 4.25 x 10~^ (32) 

From these values we can get R = 0.63, which fits the experimental data quite well. At the 
same time, the absolute branching ratios are both close to the experimental data. So we 
may conclude that when the higher twist effects are taken into account and the S-D mixing 
is considered as well, the branching ratio oi B —* iIj{3770)K can become large enough to fit 
experimental data. If a smaller value for |t| is used, the calculated decay rates are somewhat 
smaller, but still much more improved than the previous calculation. Here, the soft scattering 
effects in the spectator interactions have played an essential role. 
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Figure 2: Branching ratios of B ^ iIj{3770)K and B — > ip{3686)K (in units of 10 as 
functions of 6. The dashed hne is for B tlj{3770)K and the sohd hne for B iIj{3686)K. 

5 Discussions and Summary 

In this paper, we study the B^ tlj{3770)K~^ decay within the QCD factorization frame- 
work. If we treat '?/'(3770) as a pure l^Di state and use the leading- twist approximation for 
the kaon, we only get a very small branching ratio Bt{B ip K) =1.13 x 10~^, which is 
about 40 times lower than the experimental data. 

We further introduce the S — D mixing, combined with the calculation for the B^ 
ip{3686)K~^ decay, but still use the leading- twist approximation for the kaon, then by fitting 
the observed ratio of B^ iIj{3770)K^ to B^ iIj{3686)K^ , we find the required mixing 
angle to be about 6 = —26° or 6' = +59°. These mixing angles are not consistent with 
that obtained from other experiments. Moreover, the absolute branching ratio of B~^ — >• 
ip{3770)K~^ is still one order of magnitude lower than the experimental data. 

We then take the higher twist effects into account. By choosing proper parameters to 
characterize the end-point singularities related to the soft spectator interactions, and taking 
the S-D mixing angle to be the widely accepted value 9 = —12°, we can get a much larger 
branching ratio, and it is then possible to make the calculated rate of B~^ —>■ iIj{3770)K~^ 
close to the data. 

We would like to emphasize that in the present calculation it is the soft scattering effects 
in the spectator interactions that are essential in enhancing the decay rates, though there 
exist uncertainties for treating the soft singularities. Here, it might be useful to discuss the 
possible connection between the inclusive process B —>■ il){3770) + anything and the exclusive 
process B — > ijj{3770)K. In fact, with the nonrelativistic QCD (NRQCD) formulism |2(ij it 
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was pointed out [27] (see also j2H|) that for the D-wave charmonium inclusive production in 
B decays the color-octet cc operators in the weak decay effective Hamiltonian may play the 
dominant role by producing a color-octet cc pair at short distances, which subsequently evolve 
to a color-singlet cc (the physical charmonium) by emitting soft gluons at long distances. 
When the emitted soft gluon interacts with and is absorbed by the spectator light quark, 
the process becomes an exclusive one, such as B ^ ip{3770)K (the emitted soft gluons can 
of course hadronize into light hadrons without interactions with the spectator quark). If this 
picture makes sense, our observation in the present work that the soft scattering effects in 
the spectator interactions play the essential role in B ip{3770)K should be reasonable. 

This may also be true for the B exclusive decays involving S-wave charmonia J/iIj\22\ 
E] and ?7c|2ni, where the calculations for B J/tlj{ric)K without twist-3 soft spectator 
contributions are much smaller than the observed rates, and the enhancement effect due to 
higher twist is emphasized in^22j. For the B exclusive decays involving P-wave charmonium 
states, the situation becomes even more puzzling, that is, the measured nonfactorizable 
B XcqK decay ratejHOlISI] is large, about an order of magnitude larger than that of another 
two nonfactorizable decays B — > Xc2-^|^ and B h^K^^. These are not compatible with 
predictions based on the final state rescattering model Some of these decays are also 
studied in the PQCD approach with kt factorization ^35^, and in the light-cone sum rule 
approach|36^. In QCD factorization it is found that for the B exclusive decays involving P- 
wave charmonium states, there exist infrared divergences in the QCD vertex corrections |37| 
EH]- However, if the twist-3 soft spectator interactions dominate, we might provide a possible 
explanation for the puzzle related to B Xco{Xc2, hc)K decays, and this result will be 
presented elsewhere fSH]. 
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